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Abstract 



^ . The couphngs and decay widths of the processes (1^,2+) — > (0^,1~) + p are 

\ studied with the hght-cone QCD sum rules in the leading order of heavy quark 

CN ■ effective theory. These processes are employed to estimate the two pion transition 

. widths of the (l''',2^) doublet through the low mass tail of the p resonance. The 

\ ambiguity due to the presence of two distinct T^ states is solved. Our calculation 

00 ' shows that the two pion decay widths of the (1~'",2+) doublets are much smaller 
than the single pion decay widths. However, the Bi, B2 mesons should also have 



significant two pion decay widths around 1.5MeV. 



gn! PACS number: 12.39.Hg, 13.25.Hw, 13.25.Ft, 12.38.Lg 

. 1 Introduction 

X 

' The heavy quark effective theory (HQET) [|T| provides a systematic expansion of the heavy 

hadron spectra and transition amplitude in terms of I/tuq, where mq is the heavy quark 
mass. The spectrum pf the ground state heavy meson has been studied with the QCD 
sum rules [0] in HQET in 0]. In Q the mass of the lowest excited heavy meson doublets 
(1+, 2+) and (0"'", l"*") were studied with QCD sum rules in the heavy quark effective theory 
(HQET) up to the order of 0(l/mQ). 

The QCD sum rules was used to analyse the exclusive radiative -B-decays with the 
help of the light-cone vector meson wave function in 0. With the same formalism the 
off-shell 93* Bp and gD*Dp couplings were calculated in 0. In this work we employ the light- 
cone QCD sum rules (LCQSR) in HQET to calculate the decay widths of the processes 
(l+,2+) (0^,1") + p to the leading order of l/niq. When the LCQSR is used to 
calculate the coupling constant, the double Borel transformation is always invoked so that 
the excited states and the continuum contribution can be treated quite nicely. 

One difficult problem encountered in studying the decay widths of excited heavy mesons 
is due to the degeneracy of the spectra. There are a pair of states for any spin-parity with 
close values in their masses but quite different in magnitudes of their decay widths, except 
for the ground states. Only in the itlq 00 limit, is there a conserved quantum number 
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j£, the angular momentum of the hght component, which can be used to differentiate the 
two states. Therefore, HQET has important and unique advantage for this purpose. These 
are the motivation for our approach of using LCQSR in HQET. 

The proper interpolating current Jj'pj^^ for the states with the quantum number j, P, 
ji in HQET was given in [Q. They were proved to satisfy the following conditions 

(0|j;K'i?(0)|j',P',j;) = fpj,S,,,6pp,6^^.f--^ , (1) 

z (0|T (J5r(^)'^j'p':t'(0)) |0) = ^n'^PP'^nf^-^y ^ 9r'' ■ --g?'' 

X y dU{x - vt)Jlp^jX^) (2) 

in the mq — *■ oo limit, where rf^'"""^ is the polarization tensor for the spin j state, v is the 
velocity of the heavy quark, g^^ = g"'^ — v'^v^ is the transverse metric tensor, S denotes 
symmetrizing the indices and subtracting the trace terms separately in the sets {ai ■ ■ -aj) 
and ■ ■ ■ Pj), fpjg and Hpj^ are a constant and a function of x respectively which depend 
only on P and ji. Because of equations ([l|) and (0), the sum rule in HQET for decay 
widths derived from a correlator containing such currents receive no contribution from 
the unwanted states with the same spin-parity as the states under consideration in the 
mq oo. 



2 Sum rules for decay amplitudes 

In the present work we shall confine ourselves to the leading order of l/mg expansion. 
Denote the doublet (1^, 2+) with = 3/2 by (Si, i^g)- From covariance and conservation 
of the angular momentum of the light component in the tuq — >■ oo limit, there are three 
independent decay amplitudes for the decays of Bi, B2 to the rho meson and the ground 
states B, B* in the doublet (0^, 1^) with ji = 1/2. They are characterized by (jhj) = 
(2, 2), (1, 2), and (1, 0), where / and jh are the orbital and total angular momentum of the 
rho meson respectively P]. The S-wave amplitude is dominant. 
The amplitudes read as follows: 

M{B, ^ B*p) = I el{ie^^^P"vM^*p + ^^{^^1 ' \9L(ft)9(2,2){B^, B*) 

- VpeniqUi - \9L(ft)9(i,2){B,, B*) (3) 
+te'^"P'^r]^elv^g^i,o){Bi,B*)} , 

M{B, ^Bp)=I VpelMq^t - \qhf )9i2,^){B,, B) + gf gi,,,){B,, B)] , (4) 



M{Bl ^ B*p) = I v:ie*^{gj,,el^ - e^gijiq^qr - yrq?)9iiMBl B*) 

+r?a/3e*"e*%i,o)(52*,5*) , (5) 

M(5* ^Bp)=I r7„,„,^6'^^'^"^e;gt^.4^7(i,2)(i?2*, B) , (6) 

where r]fj_^, r]^, and e^^^ are polarization tensors for the 2"*^, l"*", 1~ heavy mesons and 
the rho meson respectively, qtfi = q^ — v ■ qVfj,. / = 1, 4^ for the charged and neutral 
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rho meson respectively. The amplitude g(2,2), 9(1,2) and 5'(i,o) in (0)-® corresponds to the 
{jhJ) = (2,2), (1,2), and (1,0) decay amplitude respectively. 

Due to the heavy quark symmetry the D-wave and S-wave amplitudes in these four 
processes are related. From vector current conservation we have another constraint. For 
example, vector current conservation requires that the tensor structure be [q^Qt — Otdt^) 
for the process Bi — > B* p. From this condition we have g(i^o){Bi, B) = — ^qfg (^2,2) {Bi, B). 
In HQET there exist only two independent coupling constants gd = g{i^2){,Bl, B) and 
9s = 9{i,o){B2, B*), corresponding to the D-wave and S-wave decays respectively. 

In order to derive the sum rules for the coupling constants g^ and gg we consider the 
correlators 

Jd^x e-^'=-(p(g)|T{j;^_ ,(0)jJ°;_3(x)}|0) = I e;{[ze^'^^'^vMQi ' HaQt) 
+ {a ^ (3)]G^,%' {00,00') + [le^^P^vMUi - HaQt) " (« f3)]G^,%\oo,oo') 

+ze^^P-v,Gf,]^; {00,00')} , (7) 



Jd^x e-'=-(p(g)|r{Jo,_,i(0)jJ_'^^ ,(x)}|0) = / e^^Kgfgf - lq^gf)Gf,]^^{oo,oo') 

+gfGf^^^^{u;,uo')], (8) 



Jd^x e-*-(p(g)|T{j;;_ ,(0)jt-|(a:)}|0) = / {{g^.^e^^ - e*,giJ{qUr 

-\grq^)G^^l!; {uo, 00') + e*-^f-Gf;7 (.., u')] , (9) 



d'x e-^'-^{p{q)\T{J,^_.{Q)jl^l7{x)m = I le^^'^^-^elqiv^qifi'^f ^^{00,00') , (10) 



^ — mi. 



where k' = k — q, 00 = 2v ■ k, 00' = 2v ■ k', ql = q^ — {q ■ v)^ and q 
The interpolationg currents are given in ^ as 

<;,3 = ^^K7'{-^){Vt - ll^ftjq , (11) 

474 = /I h^-^{^ri^r + ir^r - l9r"pt}q , (12) 

JI"-,^ = ^ h^ltQ > JI'!-,\ = \/f ^''^sg , (13) 

where is the heavy quark field in HQET and 7t^ = 7^ — v^'f). 

Let us first consider the function Gd{oo, oo') = G^^^-^ {oo, oo') in (0). As a function of two 
variables oo and oo', it has the following pole terms from double dispersion relation 

/- i/+ s.9d c c' 

^'^"^^ "^'^ = {2k_.~l'){2k^:s-uj) + 2X_.-uo' + 2A+J - a; ' ^^^^ 



where /pj^ are constants defined in (|Tp, ^^pj^ = mpj^ — mq. 
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For deriving QCD sum rules we calculate the correlator (|I0|) by the operator expansion 
on the light-cone in HQET to the leading order of l/mq. The expression with the tensor 
structure reads: 

-J^ dtj dxe-^'^-5(-a;-t;t)Tr{(z75)^(-^75)(A"-37rA)(p(g)|u(a:)rf(0)|0)} , (15) 

The rho wave function is defined as the matrix elements of nonlocal operators between 
the vacuum and rho meson state. Up to twist four the Dirac components of the rho meson 
wave function are P: 

1 







<0|^(0)a^.V^(x)|p(p,A) >=^(eWp.-e("V^)/J [ d« e'^-'^^jK /z) (16) 

(17) 



< 0|V^(0)7^V^(x)|p(p,A) >= tp^i^f^mpJ.Uu e-'^^^{<fj;{u,fi) 

+{4'^ - P,^^}fp^P fo du e-^-'^^gliu, /i) 

< 0|^(0)7m75V'(2;)|p(p,A) >= --^e^.p^e^'-^pPx'' fpTUp due-'^'i^g';{u,^i) (18) 

Due to the choice of the gauge x^Afj_{x) = 0, the path-ordered gauge factor 
Pexp {igs Jq dux^ Afj^{ux)) has been omitted. The wave functions <^J(m, yu) and ^Pp{u, jj) are 
the leading twist distributions in the fraction of total momentum carried by the quark in 
the transversely and longitudinally polarised rho meson respectively. Both wave functions 
g^jiu, jj) and g^i^u, ^) correspond to transverse spin distribution. They contain both twist- 
two and -three contributions. The twist-two contributions to the wave functions gp{u,fi) 
and gp{u, fi) can be expressed in terms of the leading twist longitudinal wave function 

j«M = i{r*^+/'A,^}, (19) 

^ 2 Jo 1 — V Ju V 

-i<,-tw,st-2(„) = 2{- r + /' d.^} . (20) 

du P ^ ' ^ Jo l-v Ju V ^ ^ ' 

All the four wave functions / = ^p^ , ip^, g^ , g"" are normalised such that /q^ f{u)du = 1. We 
shall keep only the leading twist contributions, i.e., we use the asymptotic forms for the 
following wave functions: 

(^■^(m,p) = 6m(1 -m) , (21) 

g\u,ii) = \{l + e). (22) 
g\u,^i)=Qu{l-u), (23) 
with ^ = 2m — 1. We use the wave function in [jlO| for ^'^{u) with the mock hadron mass 



replaced by the average total quark energy |[TT[] 
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^'\u) = Nexp (^- g^2^'(^_^) j [/^' + /^/^ + /^'^(l - ^)] ' (24) 

where = 0.338, = ^, ft = ^ with m = 330MeV, m = 1130MeV, and P = 320MeV. 
The wave function (E¥) satisfies the constraint from the QCD sum rule analysis of the 



moments of the rho meson wave function [T^ . The rho meson wave function is close 
to the asymptotic form, which was pointed out in p|, [1^, The rho meson couphng 
constants were determined from QCD sum rule analysis: /J = fp = 200MeV §. 
Expressing ([T5|) with the rho light-cone wave functions, we arrive at: 



G,{u:,u') = -- dt t/«e^(^-")^e»^«{/>J(«) - -/pm^l^)} . (25) 
Smilarly for the other D-wave functions in we have: 

Gf,%\uj,J) = ^G,{uj,J), (26) 

Gf,%\u,u:') = ^Ga{u,J), (27) 

Gf^%{uj,uj') = ^G,{uo,uj'), (28) 

G'^l^*{u:,uj')=Gd{uj,uj') . (29) 
For the S-wave functions, we have: 

Gf,^^;{uj,uj') = ^G,{i^,J), (30) 



Gf^^^{uj,u') = ~^G,{u,u'), (31) 



Gs{u,u') = ^G,ico,cu') = -\ml- )—^\G,{uj,uj') . (32) 
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where we define Gs((-iJ, cu') = G^^,^ [uj^uj''). 

-^G^{uj,J) = -\ml-^- J 

From the above relations we know there are two independent coupling constants only. This 
is in agreement with the discussions above. 

In order to extract we first make Wick rotation, then make double Borel transforma- 
tion to eliminate the single-pole terms in (|l^ . Subtracting the continuum contribution 
which is modeled by the dispersion integral in region lo^lo' > Uc, we arrive at: 

= le "'^ ^'^Wo{/p^^J(%)T(l - e-^) - ^/pm,(7;(«o)} , (33) 

where uq = 5^+5^' ~ T1+T2 ' -^i'-^2 are the Borel parameters corresponding to the vari- 
ables uj, uj' respectively. In obtaining (P3| ) we have used the Borel transformation formula: 
^^e"- = 5(a-i). 

For the coupling constant Qs, we have: 



^J_,i/+,3 = im^e '^T ^'^ uo{fJ^^{uo)T{l - e t ) - l/^mp^;j(Mo)} 

(/J«^J(n)T3/2(f ) - lf^m,ug^^{u)T^Mf)) , (34) 



'2 ' ' 2 
A 1 +A 3 



where /n(x) = 1 — fr is the factor used to subtract the continuum. The derivative in 

fc=0 

(0) arises from the factor {q ■ v)"^ in (|32D . We have used integration by parts to absorb the 
factor {q ■ v^. In this way we arrive at the simple form after double Borel transformation. 
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3 Determination of the parameters 



In order to obtain the coupling constants from ([33| ) and (|3J) we need to use the mass 
parameters A's and the couphng constants /'s of the corresponding interpolating currents 
as input. A„ 1/2 and /-,i/2 can be obtained from the results in [1^ as A_ 1/2 = 0.5 GeV 
and /--,i/2 — 0.25 GeV^^^ at the leading order of a^. Notice that the coupling constant 



f-,1/2 defined in the present work is a factor 1/V2 smaller than that defined in [13|. A_|_ 3/2 
is given in 0]. /+.3/2 can be determined from the formulas (34) of reference 0] derived 
from sum rules for two point correlators. The results are 

A+,3/2 = 0.82 GeV /+,3/2 = 0.19 ± 0.03 GeV^/^ , (35) 

We choose to work at the symmetric point Ti = T2 = 2T, i.e., = | as traditionally 
done in literature 0. The mass difference between (1+, 2+) and (0^, 1^) doublets is about 
0.5GeV in the leading order of HQET, which is much smaller than the large value of Ti, T2 
used below, Ti,T2 ~ 3.6GeV. So such a choice is reasonable. The values of the various 
functions appearing in (P^D, at Uq = |, are: (p'^{uo) = 1.552, gp{uo) = 1.5, (p"'^{uq) = 
-11.55 and ^"^(mq) = -12. 



4 Numerical results and discussion 

We now turn to the numerical evaluation of the sum rules for the coupling constants. The 
upper limit of T is constrained by the requirement that the continuum contribution is less 
than 40%. This corresponds to T < 2.5GeV. The lower limit of T is at the point of T 
where stability develops. This leads to T > 1.0 GeV for the sum rules (p3|) . 

Stability develops for the sum rule (g)-(|3|) in the region 1.0 GeV <T<2.5 GeV. With 
the values of pion wave functions at mq = | we obtain the sum rules for gdf^ 
Qsf- i/+ 3 as functions of T and ujc- The results are plotted as curves in Fig. 1 and 2 with 

= 3.2,' 3.0, 2.8 GeV. 

Numerically we have: 

9dL,iJ+,i = (0.17 ± 0.02)GeV2 , (36) 



Qsf-if+i = (0.095 ± 0.02)GeV^ , (37) 

where the errors refers to the variations with T and Uc in this region. The central value 
corresponds to T = 1.8GeV and Uc = 3.0GeV. 
With the values of f's in (pSj), we get: 

= (3.8±0.4±0.8)GeV-^ (38) 



= 2.1 ±0.4 ±0.4 , (39) 

where the second error takes into account the uncertainty in f's. The inherent uncertainties 
due to the method of QCD sum rules and the choice of the rho meson wave functions are 
not included here. 
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The decay width formulas in the leading order of HQET are 

r(i?' ^ Bp) = ^93I?T , 

r(i?- ^ B-p) = + 9JI9T + I^M. + + ^) g!m , (40) 



where |g| = y (mf — (m2 + sY){m\ — {m2 — s)2)/2mi, mi, m2 is the parent and decay 
heavy meson mass, s is the p meson mass and rrip = 770MeV is the p meson central mass. 

In ( ^OD the sum over charged and neutral rho meson final states has been included. In 
order to estimate the two-pion transition widths of the (l+,2+) doublet, we assume these 
transitions are dominated by the low- mass tail of the p resonance. Note 2m^ < s < mi—m2- 

We apply the leading order formulas obtained above to the excited states of charmed 
mesons assuming the HQET holds well for the charm system. The value q for these 
processes in (EOl) are calculated from the experimental mass values of the relevent particles 



Tj] . We have to smear the expression (|4^) with a Breit- Wigner form to take into account 



of the 150-MeV width of the rho meson, fis^) = --m — ^^w? — tf2- Our final results are: 

V{Di Dp) = 2.32MeV , 

T{Di D*p) = O.OSMeV , 

T{D; Dp) = O.OSMeV , 

T{D; D*p) = O.SlMeV (41) 

The masses of -81,-82 are not known experimentally. If we use tub^ = 5.755GeV and 
^B* = 5.767GeV from the quark model [§], the decay widths are: 

T{Bi Bp) = l.llMeV , 

T{Bi B*p) = 0.34MeV , 

r{B; Bp) = O.OlMeV , 

T{B; B*p) = 1.16MeV (42) 

In reference |p, the above decay widths were estimated based on the potential model 
and the assumption that the strange quark satisfies the heavy quark symmetry: 

T{Di Dp) = 6MeV , 

T{Di D*p) =< IMeV , 

T{D; Dp) =< IMeV , 

T{D; D*p) = 3MeV (43) 



r(5i ^ Bp) = 3MeV , 
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r(5i 



(44) 



In this work we have calculated the two-pion transitions of the (1+, 2+) doublets through 
the low-mass tail of the rho resonance with the method of the light cone QCD sum rules 
in the leading order of l/mg. The derived two-pion decay widths are smaller than the 
estimated values in 0]. Although we may conclude that the total width of D2 and Di is 
dominated by one pion decay, the two-pion decay channels are not negligible. Experimental 
data is not yet available for Bi and B2 mesons. From the above discussions we know that 
they should also have significant two-pion decay widths. 
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Figure Captions 

Fig. 1. The sum rule for /_ if, sgd with Uc = 3.2, 3.0, 2.8 GeV. 

'2 '^2 

Fig. 2. The sum rule for /_ if, ig^ with uj^ = 3.2, 3.0, 2.8 GeV. 
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